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PREFACE 


The  program  supported  under  this  contract  consisted  of  the 
evaluation  of  emerging  advanced  composite  materials  for  aircraft  exposed 
to  nuclear  induced  environments.  The  evaluation  focused  on  two  critical 
nuclear  response  characteristics:  (1)  thermal  flash  and  (2)  particulate 
(dust)  erosion.  Evolving  composite  materials  for  aircraft  application 
were  identified  and  various  damage  response  mechanisms  postulated. 
Materials  were  fabricated  and  various  composite  materials  were  evaluated 
analytically  and  experimentally  in  thermal  and  erosive  nulcear  generated 
environments.  The  results  of  these  analyses  are  presented  in  two 
separate  volumes  under  the  single  title  NUCLEAR  HARDNESS  SIMULATION  AND 
ANALYSIS  OF  COMPOSITE  AIRCRAFT  STRUCTURES,  Defense  Nuclear  Agency 
Contract  No.  DNA  001-82-C-0245 . 
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SECTION  1 

INTRODUCTION 

The  susceptibility  of  advanced  composite  materials  to  solid 
particle  erosion  is  becoming  an  increasing  concern.  Adler  (1981)  has 
indicated  how  laboratory  erosion  data  can  be  used  to  predict  the  amount 
of  material  removed  during  flight  through  dusty  environments.  A  very 
general  model  is  available  for  making  predictions  of  the  mean  depth  of 
penetration  (MDP)  using  relatively  simple  computational  procedures. 
This  model  is  described  in  Section  2.  Once  a  suitable  laboratory 
erosion  data  base  is  established,  estimates  of  the  MDP  for  a  wide  range 
of  flight  scenarios  and  component  geometries  can  be  made. 

BMD  Corporation  (1984)  has  also  provided  estimates  of  the  level  of 
material  removal  that  may  occur  for  advanced  composites  in  nuclear¬ 
generated  dust  environments,  however,  the  material  removal  rates  in  this 
analysis  are  based  on  results  for  ceramic  materials  and  are  not  directly 
related  to  non-metallic  composites. 

A  survey  of  the  solid  particle  erosion  literature  revealed  that 
some  experimental  data  is  available  for  polymeric  materials  and  fiber- 
reinforced,  non-metallic  composites.  Although  not  directly  applicable 
to  the  specific  composites  currently  of  interest,  the  available  data 
(summarized  in  Section  3)  is  used  to  provide  some  initial  insights  into 


the  magnitude  of  the  erosion  problem  which  may  exist  for  advanced 
composite  materials  on  aircraft  flying  through  nuclear-generated  dust. 
The  dust  distributions  have  been  modeled  by  Yoon  (1983).  Using  the 
published  erosion  data,  the  possibility  exists  for  a  significant  amount 
of  material  removal  to  occur  as  shown  in  the  examples  presented  in 
Section  4.  Additional  investigations  are  suggested  in  Section  5  to 
provide  a  more  relevant  estimate  of  the  magnitude  of  the  dust  erosion 
problem  and  the  possible  enhancement  of  the  damage  which  may  occur  by 


encounters  with  rain  environments.  The  experimental  requirements  for 
meaningful  dust  erosion  assessments  for  aircraft  are  described  in 
Section  6  where  the  capabilities  of  the  existing  facilities  are 
summarized.  A  more  general  approach  addressing  dust  erosion  effects  on 
composite  aircraft  structures  is  presented  in  Section  7.  The  magnitude 
of  the  test  requirements  and  the  analyses  which  have  to  be  accomplished 
are  clearly  indicated. 

The  lack  of  meaningful  data  for  evaluating  the  erosion  of 
composite  materials  is  a  major  obstacle  to  exercising  the  model 
described  in  Section  2.  The  estimates  provided  however  do  show  that 
there  could  be  a  significant  dust  erosion  effect  for  composite  aircraft 
components.  This  is  the  primary  conclusion  (Section  8)  resulting  from 
this  investigation. 
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SECTION  2 

EROSION  MODELING 

The  general  structure  of  an  erosion  model  which  accounts  for 
flight  through  a  particulate  environment,  the  shape  of  the  component  of 
interest,  and  the  erosion  rates  for  the  material  exposed  to  this  erosive 
environment  has  been  formulated  by  Adler  (1979a,  1981).  The  primary 
features  of  this  erosion  model  are  described.  The  general  forms  of  the 
data  required  for  implementation  of  this  erosion  model  are  clearly 
established. 

The  mean  depth  of  penetration  (MDP)  provides  a  measure  of  the 
material  removal  from  an  aircraft  component  during  flight  through  a 
nuclear  dust  environment.  The  mean  depth  of  penetration  is  determined 
as  follows: 


MDP[  t  ]  =  (1) 

pta 

where  M[ t ]  is  the  mass  removed  from  the  target  material  after  an 
exposure  time  t  in  an  erosive  environment  divided  by  the  exposed  surface 
area,  A,  and  the  density  of  the  material,  PT.  The  mass  removal  for  90 
urn  silica  sand  impacting  at  730  fps  (500  mph)  on  a  graphite  fiber 
reinforced  epoxy  target  as  a  function  of  the  exposure  time  is  shown  in 
Figure  1.  After  a  short  incubation  period  (the  nonlinear  portion  of  the 
curve)  a  steady-state  rate  of  material  removal  is  achieved.  This  is  a 
characteristic  of  the  material  removal  process  for  moderate  size 
particles  at  aircraft  flight  speeds.  The  subsequent  evaluation  of  the 
erosion  of  aircraft  components  in  nuclear-generated  dust  environments 
will  be  expressed  in  terms  of  the  steady-state  material  removal  rates. 
This  approach  provides  a  conservative  estimate  for  the  erosive  losses 
which  may  result  for  arbitrary  flight  scenarios  and  component 
geometries . 
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The  basic  measurement  of  the  MDP  defined  in  Equation  (1)  has  to  be 
generalized  for  the  flight  of  aircraft  through  airborne  dust.  The  MDP 
is  now  defined  as  the  average  depth  of  the  amount  of  material  removed 
from  the  aircraft  component  within  a  local  region  on  the  component.  The 
cumulative  effect  of  the  material  removed  from  all  of  the  local  regions 
comprising  the  entire  exposed  surface  area  of  the  component  indicates 
the  overall  shape  change  of  the  component  as  a  function  of  the  time  of 
flight  in  the  dust  environment.  This  generalized  concept  for  the  MDP 
can  be  expressed  as 


MDP  =  MDP  [X, t;  x>  rj 


(2) 


where  the  parameters  in  brackets  indicate  the  functional  dependence  of 
the  MDP.  The  coordinates  5C  denote  the  location  of  a  specific  region  on 
the  aircraft  at  a  time  t  after  the  dust  environment  is  encountered. 
The  coordinates  _x  denote  the  location  of  this  areal  element  in  space  and 
the  time  x  is  the  aging  of  the  dust  environment  at  that  location 
measured  with  respect  to  the  initiation  of  the  dust  environment.  The 
time  t  and  the  flight  path  x  can  be  selected  independently  in  order  to 
evaluate  the  level  of  erosion  damage  which  may  be  produced  by  the  range 
of  options  available  for  penetrating  the  nuclear-generated  dust 
environment . 

The  term  fm  in  Equation  (2)  is  a  representation  for  the  variations 
in  the  composition  of  the  lofted  dust  which  may  be  encountered  depending 
on  the  thermochemical  changes  which  can  occur  within  the  fireball,  the 
structural  environment  around  the  impact  site,  and  the  physics  of  the 
lofting  process.  The  influence  of  the  composition  of  the  dust  on  its 
erosiveness  is  indicated  as  a  factor  which  should  be  considered  but 
explicit  representations  for  its  effects  would  have  to  be  investigated 
in  conjunction  with  studies  of  the  dust  environment. 
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The  particle  classification  has  to  differentiate  between  angular 
or  rounded  solid  particles,  silica  sand  or  soil,  or  rain.  The 
additional  parameters  which  influence  the  rate  of  material  removal  are 
explicitly  identified  in  the  relation  which  follows. 


r\ 
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,t]  dt  (3) 


where  T  is  the  total  transit  time  and  square  brackets  are  used  to 
indicate  the  designated  quantity  is  a  function  of  the  parameters  listed 
within  those  brackets.  This  relation  simply  states  that  the  material 
removed  from  a  specific  location  on  the  exposed  surface  of  a  component 
is  equal  to  the  material  removed  per  unit  mass  of  particles  of  a 
particular  type  multiplied  by  the  total  mass  of  particles  of  that  type 
encountered  during  flight  through  a  specified  environment. 

The  conventional  measure  for  solid  particle  erosion  is  the  ratio 

of  the  mass  of  material  removed  to  the  mass  of  impacting  particles:  the 

erosion  rate.  The  erosion  rate  is  usually  determined  for  steady  state 

conditions,  that  is,  the  mass  of  material  removed  is  unchanged  for  equal 

successive  increments  of  eroding  particles.  The  erosion  rate  is 

dependent  on  the  type  of  particles  encountered,  the  particle  impact 

angle,  a  ,  the  particle  impact  velocity,  v  ,  and  the  size  of  the  impact- 
P  (k)  P 

ing  particle,  R  .  The  same  definition  of  the  erosion  rate  is 
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applicable  to  rain  erosion  except  the  erosion  rate  is  generally  a 
nonlinear  function  of  time  (or  the  number  of  drop  impacts  per  unit 
area) . 


The  designation  of  the  particle  sizes  allows  a  distribution  of 
particle  sizes  to  be  considered.  The  particle  size  distribution  curve 
can  be  divided  into  N  increments.  If  the  increments,  AR,  are  equal, 
although  this  is  not  necessary,  then  R^)  represents  the  mean  size  of 
the  range  (R  +  (k-1)  AR)  <R^)  <(R  +  kAR) .  The  effect  of  particle  size 
appears  in  the  erosion  rate  as  well  as  the  particle  concentration.  The 
cumulative  effect  on  mass  removal  is  taken  as  a  linear  superposition 
(summation)  of  the  damage  associated  with  each  particle  size.  This 
relation  can  be  modified  if  the  experimental  data  indicates  this  may  be 
required . 

A  distinction  is  made  in  Equation  (3)  between  the  vehicle's 

velocity  v  and  the  velocity  of  the  impacting  particles  v  for 

(k)  P 

particles  of  mean  radius  R  .  This  difference  acknowledges  that  the 

particle/surface  interaction  due  to  the  airflow  over  the  component  can 

modify  the  impact  velocity  of  the  particle.  The  same  considerations 

(k) 

apply  to  the  particle  impact  angle  and  the  orientation  ct  of  the 

exposed  surface. 

If  the  exposed  (local)  surface  area  is  not  planar,  it  has  to  be 
divided  into  discrete  areal  elements  due  to  the  dependence  of  mass 
removal  on  impingement  angle.  Correspondingly  Equation  (2)  can  be 
rewritten  in  terms  of  the  amount  of  material  removed  from  each  areal 
element . 

MDP  =  ^  — - — - -  MDP^1'>  (A) 

*r  A 

i 

where  MDP  ^  is  the  mean  depth  of  penetration  for  the  i^  areal  element, 
A^,  comprising  a  particular  component. 


Introducing  the  symbol  u  for  the  mass  loss  per  unit  mass  and  n  for 
the  mass  of  particles  in  a  unit  volume  (particle  concentration), 


(i)r.(i) 


MDP(±)  =  — 


.  v(lk>  R<k>]- 

pt  JtU(l)[aK)l  Lp  » 


•  n(k)  [R(k);  x ,  *t]  v(j)[x,t]t  At(j) 


The  notation  A^*^[cx^]  signifies  the  projection  of  the  area  A  ^  onto 
the  plane  normal  to  the  flight  direction.  The  time  integral  in  Equation 
(5)  is  replaced  with  a  summation.  Time  only  occurs  explicitly  in 
connection  with  the  vehicle  velocity,  since  the  steady-state  erosion 
rates  are  assumed  to  be  achieved  for  the  dust  environments  and  the 
particle  concentrations  are  a  function  of  the  time  elapsed  from  the 
onset  of  rainfall  over  a  particular  geographical  area  or  the  dust 
generated  during  a  specific  battlefield  scenario.  The  establishment  of 
the  external  environment  is  independent  of  the  vehicle's  time  of  flight. 

A  hemispherical  dome  will  be  used  to  demonstrate  how  the 
nonuniform  material  removal  from  a  nonplanar  component  can  be  obtained 
from  Equation  (3).  The  areas  in  Equation  (5)  are  evaluated  explicitly 
for  a  hemispherical  dome.  Since  the  erosion  rate  depends  on  the 
impingement  angle,  the  dome  will  be  divided  arbitrarily  into  ten  equal 
radial  segments  as  shown  in  Figure  2.  Using  the  notation  indicated  in 
Figure  2,  the  ic^  increment  of  exposed  area  of  the  target  as  a  function 
of  the  impingement  angle  A^  [a^^]  and  the  projection  of  this  areal 


normal  to  the  flight  direction  A.^  [a^ 

N 

following  equations: 


are  determined  from  the 
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A(i)  a(i)  =  2iTr2  (,cos0^l)  -  cose^) 


=  ur2  (sin29^i)  -  sin2^) 


The  numerical  values  of  Che  ten  increments  are  tabulated  in  Table  1. 

This  discrete  approach  to  the  problem  provides  a  suitable 
representation  of  the  input  data.  A  smooth  curve  can  be  drawn  through 
the  computed  MDP  data  points  as  a  function  of  dome  angle.  This  curve 
should  be  a  reasonable  approximation  to  the  actual  dome  contour  after 
exposure  to  a  defined  environment  and  flight  trajectory. 

The  flight  trajectory  is  specified  in  terms  of  the  velocity  vector 
v_.  During  the  course  of  the  vehicle's  flight  various  particle 
concentrations  will  be  encountered.  The  particle  concentrations  will  be 
functions  of  the  vehicle's  location  in  space  denoted  by  the  position 
vector  5c.  The  time-dependent  particle  concentrations  can  bo  arrayed  in 
a  three-dimensional  matrix  representing  the  space  of  general  interest  as 
determined  from  representative  dust  and  rain  distributions. 

The  position  and  velocity  vectors  for  the  vehicle  relative  to  some 
reference  coordinate  system,  such  as  the  local  cartesian  coordinates 
shown  in  Figure  3,  are 


x  =  x1e1  +  *2-2  +  *3^3 


y  =  vl£l  +  v2e,  +  v^e3 


The  components  of  the  velocity  vector  can  be  written 
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Table  1.  Geometric  parameters  associated 
with  a  spherical  dome. 


i 

r<1)/r 

B(i> 

02 

A(i)/iTr2 

4'W 

a(i> 

1 

0.1 

5°-44' 

0.010 

0.01 

87°-8' 

2 

.2 

ll°-32 ' 

.0304 

.03 

81°-22' 

3 

.3 

17°-27  * 

.0516 

.05 

75°-30' 

4 

.4 

23°-35' 

.0750 

.07 

69°-29' 

5 

.5 

30°-0' 

.1010 

.09 

63°-12 ' 

6 

.6 

36°-52 ' 

.1320 

.11 

56°-34 ' 

7 

.7 

44°-26’ 

.1720 

.13 

49°-21 ' 

8 

.8 

53°-8’ 

.2280 

.15 

41°-13 '  ! 

9 

.9 

64°-10' 

.3284 

.17 

31°-23' 

10 

1.0 

90°-0 ' 

.8716 

•19 

12°-55 ' 

11 
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=  v  cosc()  cos8 
=  v  cost})  sing 
=  v  sin<}) 


The  path  of  Che  vehicle  is  Chen  given  by 


:9+1)  -  *«>  +  v«>At 


where  che  subscripc  l  has  Che  range  1,2,3. 


The  locaCion  of  Che  vehicle  in  Che  environmenc  can  be  evaluated 


from  Equation  (11)  and  the  magnitude  of  the  velocity  vector  v^  is 
determined  from  Equation  (9).  For  numerical  evaluations  it  is 
reasonable  to  replace  the  time  integration  in  Equation  (3)  by  a 
summation  which  equals  the  total  length  of  the  path  the  vehicle  follows 
in  space  as  indicated  in  Equation  (5). 

Although  the  velocity  of  the  particles  in  a  unit  volume  at  some 
location  in  space  does  not  appear  explicitly  in  the  particle  concentra¬ 
tion  in  Equation  (5),  it  is  essential  to  know  the  average  particle 
velocity  as  a  function  of  particle  size  and  location  in  space  in  order 

C  i  k- ) 

to  evaluate  the  relative  particle  impact  speed  v  .  .  The  initial 

(k)  (k)  ^ 

particle  velocity  is  denoted  by  _v  [x.,  R  ].  The  velocity  vector  for 

(k)  ° 

the  particle  of  size  R  impacting  the  vehicle  at  a  location  x  in  space 
is  then  determined  from 


.,0)  (S)  +I(H  [*,*«]  +  vb  [v(i>,v<k),R<k)] 


where  v  is  the  velocity  of  the  vehicle  determined  from 

the  flight  trajectory, 

(k) 

v  the  velocity  of  the  particles  prior  to  interacting 

— o 

with  aerodynamic  flows  over  the  aircraft  as  determined 


from  numerical  computations  of  the  dust  distribution, 
and 


v  has  to  be  evaluated  from  an  analysis  of  the  particle 
interaction  with  the  flow  field  around  the  vehicle. 

( ilc.) 

The  impingement  angle  a  can  be  determined  from  the  velocity  vector 

( ik) 

v.  relative  to  the  surface  element  whose  normal  vector  is  in  the 

direction  tt/2  -  • 


The  expression  in  Equation  (5)  is  quite  general  and  suitable  for 
digital  analysis.  Using  a  computer  the  material  removal  rates  can  be 
evaluated  at  various  locations  on  a  component  of  irregular  shape  for 
specified  flight  trajectories  through  temporally  changing  particulate 
environments.  Once  the  vehicle  trajectory  is  established,  it  is 
necessary  to  prescribe  using  erosion  test  data  and  n  t0 

describe  the  environment.  It  is  important  to  note  that  the  erosion  data 
does  not  change;  only  the  manner  in  which  it  is  utilized  in  evaluating 
the  MDP  will  change.  Once  tests  have  been  carried  out  for  a  material  of 
interest,  for  a  range  of  impact  angles,  impact  velocities,  particle 
sizes,  and  particle  types,  this  data  base,  if  broad  enough,  can  be  used 
to  evaluate  the  erosive  effects  for  any  of  the  operational  environments 
which  may  be  encountered. 

There  are  many  factors  which  influence  the  material  removal 
process;  however  the  particle  size  exponent  and  the  velocity  exponent 
will  be  of  primary  concern,  since  these  parameters  are  common  to  the 
majority  of  the  particle  erosion  analyses  which  can  be  found  in  the 
literature.  The  lack  of  agreement  between  the  results  of  the  theories 
based  on  kinetic  energy  exchange  between  the  impacting  particle  and  the 
target  has  been  a  perplexing  issue  in  the  field  of  solid  particle 
erosion.  The  general  form  of  the  expression  for  volume  removal  (volume 
removed  per  particle  impact)  can  be  written, 

=  g  v3  Rbf  (  ;t  )  (13) 
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where  g  is  a  functional  representing  the  dependence  on  all  of  the 
parameters  not  explicitly  stated  and 

v  is  the  particle  impact  velocity 

R  is  the  mean  radius  of  the  impacting  particle 

f(ct  )  is  a  function  of  the  attack  angle  a 
p  P 

a  is  the  velocity  exponent 

b  is  the  particle  size  exponent. 

The  velocity  exponent  has  been  generally  found  to  range  from  two 
to  three  for  materials  displaying  ductile  behavior  for  a  wide  range  of 
particle  impact  conditions,  however  exponents  (up  to  4)  have  been 
obtained  by  Grant  and  Tabakoff  (1975)  and  Tabakoff,  et  al  (1979). 
Tabakoff 's  results  are  based  on  detailed  observations  of  the  aerody¬ 
namics  of  the  particle  flow  and  the  rebound  characteristics  of  the 
particles  impacting  the  eroding  surface  of  the  target.  The  velocity 
exponents  for  brittle  materials  range  from  1  to  6.5  (cp.  Gulden,  1979a; 
Sargent,  et  al.,  1979).  The  velocity  exponent  depends  on  the  size, 
material,  shape,  and  angle  of  attack  of  the  impacting  particle,  however 
it  does  not  appear  to  be  significantly  affected  by  the  microstructure  of 
the  target  material. 

The  particle  size  exponent  has  not  been  adequately  evaluated.  It 

is  generally  stated  that  the  volume  removal  in  metals  is  independent  of 

particle  size  for  particle  sizes  in  excess  of  100  Urn  (Tilly,  1969).  The 

variations  in  particle  size  effects  have  already  been  noted  for  smaller 

particle  sizes  in  brittle  materials  (Sheldon  and  Finnie,  1966). 

Gulden's  test  data  (Gulden,  1979a, b)  shows  b:r4  for  magnesium  fluoride 

impacted  by  quartz  particles,  however  hot-pressed  silicon  nitride 

impacted  by  silicon  carbide  particulates  displayed  a  b=4  dependence 

while  impacts  by  quartz  produced  a  b*3  dependence  as  well  as  a  change  in 

a=*4  to  a=l  for  particle  sizes  from  8  to  940  Um.  A  definite  threshold 

for  achieving  an  V  ^  relation  was  found  for  reaction-bonded  silicon 

o 

nitride  impacted  by  quartz  particles.  The  above  results  are  for  normal 


imoacts  on  the  target  surface.  The  dependence  of  the  particle  size 
exponent  on  the  particle  material  is  evident,  but  it  is  not  known  how  it 
would  vary  with  the  angle  of  attack. 

The  complete  characterization  of  a  particular  impact  condition 
(particle  impact  type  and  target  material)  involving  the  particle  size, 
impact  velocity,  and  impingement  angle  requires  numerous  erosion  tests. 
Thus  a  reliable  experimental  data  base  for  the  development  as  well  as 
comparison  of  analytical  results  is  rarely  found  in  the  particle  erosion 
literature  which  leaves  many  aspects  of  the  subject  unsettled.  The 
analytical  approaches  which  have  been  undertaken  to  describe  material 
removal  due  to  both  solid  and  liquid  particle  erosion  have  been 
described  and  evaluated  by  Adler  (1979b, c). 


SECTION  3 


SOLID  PARTICLE  EROSION  OF  COMPOSITE  MATERIALS 

Experimental  evaluations  of  the  material  removal  rates  for 
polymeric  and  composite  materials  are  summarized  in  Table  2.  This  is 
all  of  the  data  which  could  be  found  in  the  open  literature.  Some 
general  observations  concerning  this  data  are  provided. 

A  somewhat  more  extensive  data  base  exists  for  rain  erosion 
evaluations  of  polymeric  materials  than  for  solid  particle  erosion. 
However  it  is  important  to  note  that,  while  attempts  have  been  made  to 
correlate  the  data  generated  in  each  of  these  environments,  little 
correspondence  exists  between  the  respective  material  removal  rates. 
Thus,  based  on  the  author's  experience  related  to  both  solid  particle 
and  rain  erosion  evaluations,  it  is  not  possible  to  use  rain  erosion 
data  to  predict  the  erosion  rates  for  solid  particle  impacts  (Adler, 
1979a). 

Setting  this  issue  aside,  the  solid  particle  erosion  investiga¬ 
tions  summarized  in  Table  2  will  be  used  as  the  data  base  for  estimating 
the  material  removal  rates  for  potential  nuclear  dust  environments.  The 
recent-  work  of  Zahavi  and  Schmitt  ( 1981a, b)  is  of  interest,  but  the 
particle  impact  velocities  (about  42  ms  )  are  too  low  to  provide 
relevant  material  removal  rates. 

There  are  some  general  trends  in  the  solid  particle  erosion  data 
which  remain  to  be  established  for  fiber-reinforced  composites.  For 
ductile  metals  the  steady-state  erosion  rates  become  constant  with 
respect  to  particle  size  beyond  a  certain  particle  size  (on  the  order  of 
100  am);  the  critical  particle  dimension  is  dependent  on  the  particle 
type  and  impact  velocity  (Goodwin,  Sage  and  Tilly,  1969;  Tilly  and  Sage, 
1970).  This  limitation  does  not  appear  to  be  applicable  to  glass  which 
is  representative  of  brittle  response.  The  maximum  (steady-state) 
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Table  2.  Summary  of  test  data  for  solid  particle  erosion 
of  composite  materials. 


Reference 


(1974) 


Tilly  and  Sage 
(1970) 


Tilly  69) 


.er  1*73) 


Car.avi  and 
Schmitt  f 1981a* 


-ahavi  and 
Schmitt  ‘’1981b) 


Particle 

Particle 

Test 

Velocitv 

Impact 

Material 

Material 

SUe 

Apparatus 

Range 

Angle 

(mm) 

(ms-1) 

CFRP  180/377. 

Si02- 

0.2  to  0.25 

Rotating  arm 

100-400 

30  to  90 

Astrocoat  on  GFRP 

Si02- 

0.1  to  0.125 

90° 

0.2  to  0.25 

250 

30  to  90 

0.3  to  0.4 

90“ 

GFRP  +  400  uo 

SiO  - 

0.2  to  0.25 

250 

90“ 

Fluorocarbon  coating 

GFR  epoxy 

SiO.-rounded 

0.044  to 

unknown 

47-109 

90“ 

(unidirectional 

0.250 

graphite) 

SiO^-rounded 

0.077  to 

70 

30  to  90 c 

0.125 

Polyurethane 

SiO^-angular 

0.125  to 

Rotating  arm 

123 

90° 

0.150 

N'ylon  (Type  66) 

60-305 

90“ 

Glass-reinforced  nylon 
30%  glass  fiber  in 
Type  66 

Carbon-reinforced  nylon 
25%  carbon  fiber  in 
Type  66 

Fiberglass 

70%  glass  in  epoxy 

Polypropylene 
.Vvlon  (Type  66) 

Carbon-reinforced  nylon 
25%  carbon  fiber  in 
Type  66 

Class-reinforced  nylon 

Fiberglass 
7 OX  glass  in  epoxy 

Polypropylene 

;fr 

U  v'o 

Epoxy  "unidirectional) 


Si0„ -angular  0,06  Co 

0.125 


Air  blast 


I  SiGn -angular 


Hass-reinforced  polv-  SiO., -rounded  0.210  to 


Rotating  am 


Air  blast 


-elongated 


Class- re info  reed 
polyimide 

Class-epoxv  laminates 
Typically  35%  resin, 
65%  reinforcement 

Coatings  - 
Polyurethane 
■'MIl-C-83236) 

Polyurethane 
' MIL-C-33231) 

Fluorocarbon 

' AF-C-VSW-15-15) 


0.210  to 

0.297 


Air  blast 


40°,  90°, 
70  to  90°  ' 


30°  to  90° 


15 4  to  90 3 
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erosion  race  for  ductile  metals  occurs  at  an  impact  angle  around  20°  , 
while  normal  (90°)  particle  impacts  are  the  most  damaging  for  brittle 
materials.  For  very  small  particles  (c  10  Pm)  the  maximum  erosion  rate 
for  brittle  ceramics  shifts  to  impact  angles  less  than  90°  (Sheldon  and 
Finnie,  1966). 

The  trends  in  the  erosion  rates  for  polymeric  composites  are 
dependent  on  the  specific  composite  material  being  evaluated  and  are 
typically  intermediate  to  the  general  trends  observed  for  ductile  and 
brittle  materials.  For  polymeric  composites  there  is  a  tendancy  for  a 
mass  increase  for  certain  impact  conditions  (as  normal  incidence  is 
approached)  due  to  particle  embedding  rather  than  mass  removal  during 
the  initial  stages  of  the  erosion  process. 

Details  of  the  experimental  apparatus  and  procedures  used  are 
relatively  scant  in  the  published  reports  along  with  the  test  procedures 
used.  The  limited  information  that  is  available  will  be  used  to 
reconstruct  the  manner  in  which  the  stated  results  were  obtained. 

Tilly  and  his  co-workers  (Tilly,  1969;  Tilly  and  Sage,  1970)  used 
both  a  sand  blast  apparatus  and  a  rotating  arm  test  configuration  as 
indicated  in  Table  2.  The  rotating  arm  provided  impact  velocities  more 
representative  of  those  encountered  in  gas  turbines  operating  in  dusty 
environments,  while  the  sand  blast  apparatus  was  used  to  obtain  elevated 
temperature  data.  A  single  specimen  was  mounted  at  each  end  of  a  10.16 
cm  (4  in.)  arm  and  for  higher  velocities  (approaching  550  ms  ^  )  a  15.24 
cm  (6  in.)  arm.  Sand  was  dropped  from  a  dispenser  at  one  point  along 
the  circumference  of  the  circle  through  which  the  specimen  passed.  The 
height  of  fall  of  the  sand  could  be  adjusted  so  that  all  of  the  sand 
impacted  the  specimens.  Tilly  indicates  that  at  most  7%  of  the  sand 
released  may  miss  the  specimens.  The  specimens  were  1  cm  (0.4  in.) 
squares  machined  with  a  shank  for  attachment  to  the  rotating  arm.  It  is 
assumed  that  the  polymeric  samples  were  fastened  to  a  metal  backing  in 
order  to  obtain  mass  loss  data  for  these  materials.  In  the  sand  blast 


device  the  particle-laden  air  stream  was  directed  at  the  central  portion 
of  the  specimen.  The  quantity  of  particles  impinging  the  specimen  is 
taken  to  be  proportional  to  the  projected  area  of  the  test  specimen 
normal  to  the  flow. 

The  AFML/Bell  rotating  arm  erosion  facility  (Adler,  1973)  has  a 
rotating  arm  which  places  the  specimen  at  a  radial  distance  slightly 
greater  than  2.75  m  (9  ft).  The  short  arm  radius  used  by  Tilly  did  not 
let  the  particles  fall  the  width  of  the  specimen  before  they  were 
impacted  again.  The  large  radius  of  the  Bell  arm  allows  the  particles 
to  fall  more  than  30  cm  before  they  are  impacted  on  the  next  revolution 
of  the  arm.  The  Dornier  System's  rotating  arm  (Behrendt,  1974)  has  a 
radius  of  1.2  m  (3.94  ft).  No  information  is  provided  in  the  paper  by 
Williams  and  Lau  (1974)  as  to  the  apparatus  or  procedures  which  were 
used  to  obtain  the  data  reported. 

The  actual  number  and  condition  of  the  particles  striking  the 
specimens  in  these  erosion  facilities  are  typically  determined  indi¬ 
rectly.  Adler  (1973)  calculated  the  number  of  impacts  for  the  AFML/Bell 
rotating  arm  using  the  following  relation  for  the  test  arrangement  shown 
in  Figure  4, 

N  =  - ~= ;  yr  d’  COS0  (14) 

720  V  2gh 


where  is  the  number  of  particle  impacts  the  specimen  experiences  per 
revolution  of  the  rotating  arm, 

N  is  the  number  of  particles  released  into  the  path  of  the 
specimen  per  minute, 

\/ 2gh  =  21  .233  fps  for  the  AFML/Bell  erosion  facility, 
l  is  the  lateral  dimension  of  the  sand  ejection  nozzle, 


V  is  the  lateral  dimension  of  the  exposed  area  on  the  specimen's 
surface  in  inches. 


d  is  the  transverse  dimensions  of  the  exposed  area  on  the 
specimen's  surface  in  inches. 

The  specimen  used  in  Adler's  experiment  was  a  half  cylinder:  1.27  cm 
(0.5  in.)  diameter  x  5.08  cm  (3  in.)  long.  The  erosion  rates  with  this 
geometry  represent  a  range  of  impact  angles  which  distorts  the  values 
obtained.  For  normal  impacts,  9  =  0°  as  defined  here,  l  »  1.375  in.,  V 
=  2.5  in.,  and  d'  =■  0.5  in.  Since  £'  is  greater  than  l,  the  ratio  i/l' 
is  taken  to  be  unity;  however,  the  particle  spray  expands  as  it  exits 
from  from  the  nozzle,  so  deviations  from  the  ideal  conditions  do  occur. 
The  nominal  90  um  quartz  sand  was  dispensed  at  a  rate  of  260  g/min.  The 
particle  size  distribution  is  shown  in  Figure  5  and  is  skewed  to  the 
larger  particle  dimensions.  Assuming  an  equivalent  spherical  diameter 
of  90  ttm  and  a  density  for  quartz  of  2.2  g/cc,  the  value  of  N  in 

O 

Equation  (12)  is  determined  to  be  3.10  x  10  particles/min.  Using  these 

values  in  Equation  (12)  yields  a  value  of  =  9.30  x  10^  particles/rev. 

At  v  =  222  ms  1  (730  fps)  the  RPM  for  the  AFML/Bell  arm  is  1.05  v  or 
P  P 

7t>6.  The  particle  impacts  per  second  are  then  estimated  to  be  1.19  x 
10^.  This  value  can  be  used  with  the  measured  mass  loss  per  exposure 


increment  in  Figure  1  to  obtain  an  approximation  of  the  erosion  rate. 
It  is  again  important  to  note  that  the  specimen  was  a  half-cylinder  and 
so  the  erosion  rate  in  Table  2  cannot  be  reliably  associated  with  a 


single  impingement  angle.  With  this  reservation 
assigned  for  the  initial  evaluation. 


=  90°  will  be 


Behrendt  (1974)  does  not  provide  explicit  information  on  how  the 
mass  of  impacting  particles  was  obtained,  so  it  is  not  possible  to 
evaluate  the  consistency  in  the  erosion  rates  he  obtained  and  those 
obtained  by  other  investigators  using  rotating  arms.  Behrendt 's  erosion 
evaluations  are  moderately  extensive  providing  the  variations  of  erosion 
rates  with  impact  velocity  and  impact  angle  as  shown  in  Figure  6  to  8 
for  a  graphite  fiber  reinforced  composite  (the  exact  constituents  are 
not  specified).  The  samples  are  cylindrical  with  a  16.8  mm  (0.66  in.) 
diameter  and  a  thickness  of  2.5  to  5.0  mm  (0.1  to  0.2  in.). 
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Figure  5.  Particle  size  distributions  for  70-micron  glass 
beads  and  90-micron  sand  (Adler,  1973). 
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The  dependence  of  the  erosion  rate  on  the  impact  parameters  was 
defined  in  Equation  (13). 

U  =  g  va  Rb  f(ct  ) 

6  p  P 

where  g  is  a  functional  representing  the  dependence  on  all  of  the 
parameters  not  explicitly  stated,  a  is  the  velocity  exponent,  and  b  is 
the  particle  size  exponent.  Behrendt's  results  indicate  that  a  s=2.4 
and  the  functional  dependence  on  at  250  ms  ^  can  be  obtained  from 
Figure  8.  The  maximum  erosion  rate  occurs  at  an  impact  angle  of  45°. 
On  the  other  hand,  Williams  and  Lau  (1974)  find  from  their  experiments 
that 

’0  =  0.01  (2R)°*9  [v  sin  (ap  +  24°)  ]3’4  (15) 

This  expression  is  based  on  100  Mm  quartz  particles  impacting  unidirec¬ 
tional  graphite  fiber  reinforced  epoxy  speciment  at  70  ms  1  (230  fps). 
The  maximum  erosion  rate  for  the  test  conditions  used  by  Williams  and 
Lau  occurs  at  an  impact  angle  of  65°  (Figure  9).  It  is  also  important 
to  note  that  the  velocity  exponent  in  Equation  (15)  is  3.4  which  is 
considerably  higher  than  that  found  by  Behrendt  (1974)  but  the  impact 
velocity  is  quite  low. 

Tilly  and  his  co-workers  (Table  2)  find  that  the  erosion  rate  for 
nylon  and  glass  carbon  fiber  reinforced  nylon  has  a  maximum  at  an 
impingement  angle  of  30°  (Figure  10).  Quartz  particles  ranging  from  60 
to  125  Mm  sieve  diameters  imparted  the  specimens  at  v^  =  104  ms  (341 
fps).  Tilly  demonstrated  that  a  significant  mass  of  quartz  particles 
had  to  impact  the  polymeric  test  samples  before  the  erosion  rates  would 
begin  to  approach  their  steady-state  values.  This  effect  is  greater  for 
90°  impacts  than  for  40°  collisions  (Figure  11).  400g  of  quartz 

particles  were  used  to  obtain  the  erosion  rates  in  Figure  10,  however 
according  to  Tilly's  evaluation  of  the  influence  of  the  mass  of  the 
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Figure  10.  Variation  of  erosion  rates  for  polymeric  materials 
with  impact  angle  (Tilly,  1969). 


a)  90°  Impacts  b)  40°  Impacts 


Figure  11.  Variation  of  weight  loss  with  quantity  of  abrasive 
for  polypropylene  and  nylon  (Tilly,  1969). 


impacting  particles,  the  erosive  response  is  still  not  stabilized  for 
400g  of  quartz  particles.  This  observation  raises  doubt  as  to  the 
accuracy  of  the  data  reported  in  Figure  10.  In  addition  the  data  points 
are  not  sufficient  to  confirm  that  the  results  for  nylon,  carbon- 
reinforced  nylon,  and  glass-reinforced  nylon  are  nearly  identical  in  the 
vicinity  of  a  =  30°  (Figure  10).  For  a  fairly  broad  range  of  materials, 
Tilly  and  Sage  (1970)  found  that  the  velocity  exponent  was  approximately 
2.3  for  normal  impacts  of  125  to  150  hm  quartz  particles  impacting  from 
60  to  300  ms  ^  as  can  be  seen  in  Figure  12. 

The  erosion  rates,  ’J  ,  which  can  be  determined  from  the  publica¬ 
tions  described  above  are  summarized  in  Table  3  for  the  particle  impact 
conditions  which  are  of  interest  for  assessing  aircraft  vulnerability  to 
nuclear-generated  dust.  Although  a  sizable  increase  is  seen  in  the 
erosion  rates  at  intermediate  impact  angles  at  low  impact  velocities  (as 
shown  in  Figure  10) ,  the  difference  between  the  maximum  erosion  rate  at 
an  intermediate  angle  compared  to  normal  collisions  may  not  be  as 
drastic  at  impact  velocities  in  the  vicinity  of  Mach  1  (as  indicated  in 
the  plot  in  Figure  8).  The  values  of  u  for  normal  impacts  in  Table  3 
will  therefore  be  assumed  to  be  slightly  low  but  indicative  of  the 
maximum  erosion  rates. 

In  order  to  convert  the  volumetric  erosion  rates  reported  by  Tilly 
and  Sage  (  1970)  to  mass  loss  per  unit  mass  of  impacting  particles,  the 
densities  of  the  composites  they  evaluated  have  been  estimated.  These 
estimated  densities  are  indicated  in  Table  3.  The  composite  materials 
are  all  quite  different,  so  the  available  data  can  only  be  used  to 
provide  a  very  crude  estimate  of  the  solid  particle  erosion  rates 
current  composite  materials  may  experience  for  better  defined  erosive 
conditions.  The  erosion  rates  from  Tilly  and  Sage  (1970)  show  the  range 
of  values  that  may  occur,  however  the  graphite-fiber  reinforced  epoxy 
data  obtained  by  Adler  (1973)  and  Behrendt  (1974),  which  may  be  more 
representative  of  present  day  composites,  will  be  used  to  obtain  some 
idea  of  the  MDPs  for  aircraft  flying  through  nuclear-generated  dust 
environments . 
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Figure  12.  Velocity  dependence  of  composites  tested  with 

125-150  um  quartz  at  90°  (Tilly  and  Sage,  1970) 
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Table  3.  Erosion  rates  for  quartz  particles  impacting 
composite  materials  at  normal  incidence. 


Source 

Material 

Particle 

Impact 

Size 

Velocity 

(mm) 

(ms--*-) 

1 

Tilly  and  Sage 

30%  glass  fibers  in 

125-150 

185 

4.93 

(1970) 

type  6,6  nylon 

245 

10.9 

p  ~  1.56  g/cc. 

305 

17.14 

70%  glass  fibers  in 

185 

53.5 

epoxy  (fiberglass) 
p  s  2.05  g/cc. 

245 

109. 

305 

170. 

25%  carbon  fibers  in 
type  6,6  nylon 

185 

245 

4.05 

9.00 

p  a  1.28  g/cc 

305 

14.06 

Type  6,6  nylon 

185 

1.28 

p  s  1.16  g/cc 

245 

2.67 

305 

4.53 

Behrendt  (1974) 

GFRP  180/37% 

200-250 

100 

6.08 

200 

26.4 

250 

43.5 

300 

70.8 

400 

137.7 

Adler  (1973) 

31  v/o  graphite 
in  epoxy 

90 

222 

19.10 

The  dust  distributions  (Yoon,  1983)  are  specified  in  terms  of  the 
particle  size.  Unfortunately  the  erosion  rates  obtained  by  Adler  (1973) 
and  Behrendt  (1974)  are  for  a  single  mean  particle  size,  so  it  is  quite 
difficult  to  specify  the  variations  of  u  with  particle  size.  Several 
forms  of  the  dependence  of  U  on  particle  size  are  likely  to  exist  based 
on  the  trends  in  the  solid  particle  erosion  response  of  other  materials. 
The  available  data  is  plotted  in  Figure  13  as  a  function  of  the  mean 
particle  size.  The  range  in  the  form  of  the  dependence  of  t)  on  particle 
size  is  indicated  in  Figure  13.  At  this  time  it  is  not  possible  to 
select  one  form  over  another,  but  the  data  deficiency  could  be 
eliminated  with  a  small  experimental  effort.  It  was  stated  at  the 
beginning  of  this  section  that  the  steady-state  erosion  rates  for 
ductile  metals  become  constant  with  respect  to  particle  size,  however 
the  erosion  rates  for  more  brittle  materials  may  not  display  this  same 
effect  but  may  tend  to  increase  for  a  considerably  larger  range  of 
particle  dimensions.  These  two  possibilities  are  reflected  in  the 
erosion  rates  for  the  smaller  particle  dimensions. 

The  solid  particle  erosion  rates  for  non-metallic  composites 
obtained  from  the  existing  literature,  be  that  what  they  may,  will  be 
used  in  Section  4  to  provide  very  preliminary  estimates  of  the  removal 
rates  for  hypothetical  aircraft  trajectories  traversing  the  continental 
United  States. 


SECTION  4 


MASS  REMOVAL  ESTIMATES 

The  extent  of  the  damage  which  may  result  from  flight  through 
nuclear-generated  dust  will  be  estimated  using  the  erosion  model 
described  in  Section  2  and  the  extrapolations  of  the  experimental 
erosion  rates  summarized  in  Section  3.  The  specific  conditions  to  be 
evaluated  are  the  hypothetical  deployment  of  aircraft  from  the 
southwestern  United  States  at  various  times  after  a  nuclear  attack  on 
missile  sites  within  the  continental  United  States. 

The  general  relation  in  Equation  (5)  will  be  evaluated  using  the 
dust  distributions  calculated  by  Yoon  (1983).  Representative  dust 
distribution  plots  are  available  as  shown  in  Figures  14  to  16.  Two 
trajectories  passing  through  the  dust  fields  have  been  selected  for 
aircraft  flying  across  the  continental  United  States  without  consider¬ 
ation  of  how  realistic  they  might  be.  These  trajectories  will  be  used 
to  demonstrate  how  the  temporal  development  of  the  dust  fields  can  be 
introduced  into  the  erosion  model  to  evaluate  the  MDP.  The  two 
trajectories  labeled  A  and  B  are  superposed  on  the  dust  distributions  in 
Figures  14  and  16.  Since  only  a  limited  number  of  time  increments  have 
been  computed  for  the  dust  distribution,  the  evaluation  of  Equation  (3) 
has  to  be  carried  out  for  a  much  larger  time  increment  (essentially  a 
static  dust  distribution)  than  would  be  the  case  if  the  model  could  be 
adequately  implemented.  A  more  reasonable  time  increment  for  these  dust 
distributions  would  be  At  *  100  sec.  for  an  aircraft  flying  at  300  to 
400  ms  ^  .  The  aircraft  would  fly  roughly  20  nmi  during  this  time 
increment  for  a  flight  speed  of  375  ms  ^  (c.750  knots).  It  would  also 

be  desirable  to  have  the  dust  concentrations  isocontours  at  intervals  of 

3 

3,  15,  30,  45,  60,  90,  120,  and  greater  than  150  mg/m  . 

The  dust  field  calculations  (Yoon,  1983)  provide  the  mean  and 
range  of  particle  sizes  in  the  atmospnere  as  a  function  of  time  after 


laydown  as  indicated  in  Figure  17.  For  the  restricted  conditions  which 
can  be  investigated  at  the  present  time.  Equation  (5)  reduces  to 


=  -f-  y;u-n(j)-L(j) 

0 


Normal  impacts  are  assumed  so  there  is  no  angular  dependence.  The 
erosion  rate  is  then  simply  a  function  of  the  particle  size  and  impact 
velocity.  The  particle  size  is  determined  from  Figure  17  for  the  values 
of  time  after  laydown  and  therefore  is  assigned  a  fixed  value. 
Correspondingly  the  mean  erosion  rate  remains  unchanged  which  would  not 
be  the  case  if  the  time  interval  for  the  dust  distribution  was  reduced 
to  15  min.  For  lack  of  a  flight  plan  the  aircraft  velocity  is  also  held 


constant,  so  from  Equation  (5)  J  vUJ  j  At^j;  =»  L'' 
aircraft  travels  through  dust  with  concentration  n( 
the  generic  composite  is  assumed  to  be  1.65  g/cc. 


the  distance  the 
.  The  density  of 


The  MDP  for  trajectories  A'  and  B  can  now  be  evaluated  for  the 
conditions  stated  above.  The  MDPs  for  a  range  of  conditions  are  listed 
in  Table  4  for  trajectory  A  and  in  Table  5  for  trajectory  B.  Referring 
to  Tables  4  and  5,  it  is  seen  that  the  amount  of  material  removed  for 
both  the  short-  and  long-time  distributions  is  on  the  order  of  1  mm  (40 
mil).  The  MDP  values  can  vary  widely  from  this  value  due  to  the 
uncertainties  in  the  erosion  data  as  well  as  the  dust  distribution 
calculations.  Also  no  effort  was  made  to  make  worst  case  predictions. 
Flight  through  higher  dust  concentration  zones  would  increase  the 
removal  rates  several  fold. 
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Table  4.  Mean  depth  of  penetration  estimates  for 
flight  trajectory  A. 


SECTION  5 


DISCUSSION 

The  relevant  erosion  daca,  five  experimental  data  points  in  Figure 
13,  for  two  different,  ill-defined  graphite-fiber  reinforced  composites 
has  been  used  as  the  basis  for  constructing  highly  speculative  rela¬ 
tions  between  the  erosion  rates  and  the  particle  dimensions.  This  data 
is  required  since  the  mean  particle  dimensions  in  the  dust  environment 
predictions  are  a  function  of  the  time  after  cloud  formation  (Figure  j 

17).  The  erosion  rates  recorded  in  Figure  13  can  vary  considerably  from 
the  stated  values  due  to  the  high  degree  of  uncertainty  in  the  test 
procedure  and  the  quality  of  the  test  specimens  compared  to  current 
composite  materials.  However  the  information  provided  is  indicative  of  j 

the  prevailing  state-of-the-art  for  solid  particle  erosion  of 

nonmetallic  composites  and  coated  materials  for  impact  velocities  in  the 
vicinity  of  Mach  1.  Figure  13  represents  the  best  estimates  that  can  be 
made  without  initiating  additional  testing  to  remove  the  prevailing  data 
gap.  The  dust  environment  predictions  also  contain  a  high  level  of 
inherent  uncertainty  which  could  significantly  influence  the  magnitude 
of  the  estimated  material  removal  rates.  In  view  of  all  of  these 

factors  and  additional  considerations  which  will  be  discussed  shortly, 

I 

it  is  evident  the  detrimental  effects  of  the  dusty  environments  which  1 

may  be  encountered  by  composite  aircraft  components  cannot  be  ignored. 

The  90°  particle  impact  data  used  to  make  the  MDP  estimates  in 
Tables  4  and  5  may  not  be  the  worse  case  condition.  Based  on  the  i 

experimental  erosion  rates  for  various  composite  materials,  it  appears 
that  the  critical  impact  angle  (the  angle  at  which  the  maximum  erosion 

rate  occurs)  is  dependent  on  the  type  of  composite  material  being 

evaluated.  The  critical  impact  angle  can  range  from  30°  to  90°  (Figure 
10).  Maximum  erosion  rates  are  seen  to  occur  at  45°  (Figure  8)  and  65° 

(Figure  9).  Thus,  until  the  erosion  rates  are  measured  for  a  specific 
composite  construction,  there  is  no  assurance  that  non-normal 


impingement  angles  will  be  beneficial  in  alleviating  the  severity  of  the 
amount  of  material  removed. 

The  material  removal  evaluated  in  Tables  4  and  5  is  for  aircraft 
velocities  in  the  vicinity  of  Mach  1;  however  if  the  flight  velocity 
increases  to  Mach  2,  the  material  removal  is  increased  by  a  factor  of 
5.28  according  to  the  power  law  relation  for  the  erosion  rates  obtained 
from  Behrendt's  (1974)  data.  Extrapolation  to  the  higher  velocity 
assumes  that  no  significant  change  occurs  in  the  mode  of  material 
removal  between  these  two  velocities. 

The  erosion  data  for  the  graphite-reinforced  polymers  listed  in 
Table  3  was  used  to  provide  the  estimates  of  mass  removal  in  Tables  4 
and  5.  If  the  erosion  data  for  a  material  with  the  characteristics  of 
fiberglass  was  used  instead,  the  mass  removal  could  conceivably  increase 
by  a  factor  of  2.5  for  the  normal  impact  condition.  On  the  other  hand, 
a  more  resilient  matrix  material  than  epoxy  (the  nylon  composites  in 
Table  3)  could  help  to  reduce  the  mass  removal.  The  variation  of  more 
than  an  order  of  magnitude  in  the  experimental  erosion  rates  for 
different  composite  materials  tends  to  reinforce  the  need  for  estab¬ 
lishing  reliable  erosion  data  for  a  specific  composite  construction. 

The  possibility  exists  for  aircraft  which  have  been  exposed  to 
dust  to  subsequently  encounter  rain  along  their  flight  path.  Although 
it  may  exist,  no  experimental  data  could  be  found  for  providing  insights 
into  the  synergistic  effects  of  the  sequential  imposition  of  a  dust  and 
rain  environment  on  any  material.  The  interaction  between  the  erosion 
characteristics  for  each  of  these  particle  types  could  be  quite  strong 
and  warrants  consideration  if  the  potential  for  damage  is  to  be 
representative  of  a  broad  range  of  flight  scenarios.  While  no  directly 
applicable  erosion  data  is  available,  conjectures  can  be  made  as  to  the 
extent  of  the  material  loss  based  on  a  knowledge  of  the  erosion  behavior 
of  materials  individually  exposed  to  these  two  erosive  environments.  In 
particular,  Adler  (1973)  evaluated  the  rain  erosion  behavior  of  the 
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graphite-reinforced  epoxy  material  in  Table  3  as  well  as  its  solid 
particle  erosion  characteristics.  The  rain  erosion  data  is  shown  in 
Figure  18.  The  rain  field  is  calibrated  to  1 .8  mm  mean  diameter  water 
drops  falling  at  the  rate  of  2.54  cm/hr.*  It  is  seen  that  some  time 
(the  incubation  period)  is  required  before  any  rapid  rate  of  material 
removal  takes  place.  Finally  a  steady-state  rate  of  removal  is 
achieved.  Using  the  information  in  Adler's  report,  the  number  of  drops 
impacting  the  specimen  and  correspondingly  the  mass  of  water  impacting 
the  surface  can  be  approximated.  The  steady-state  erosion  rate  (for  1.8 
mm  diameter  drops)  was  evaluated  from  Figure  18  and  was  found  to  be  52 
mg/g  compared  to  19.1  mg/g  for  the  solid  particle  environment.  The  same 
specimen  geometry,  a  half  cylinder,  was  used  in  the  rain  erosion  tests 
also.  If  the  region  of  the  leading  edge  in  which  the  drop  impacts 
normal  to  the  surface  of  the  specimen  was  used  in  the  above  erosion  rate 
evaluation  instead  of  the  projected  area  of  the  specimen,  the  erosion 
rates  would  be  significantly  higher  (possibly  by  as  much  as  a  factor  of 
four) . 

The  general  character  of  the  damage  is  distinctly  different  in 
each  case.  The  multiple  solid  particle  impacts  produce  a  relatively 
smooth  eroded  surface  with  each  particle  removing  a  small,  localized 
volume;  the  water  drops  can  penetrate  the  bulk  of  the  specimen  once 
cracks  and  pits  are  formed  on  the  surface  and  remove  relatively  large 
volumes  of  material  by  a  surface  penetration  and  upheaval  mechanism 
(Adler,  1979b). 

Returning  to  the  potential  for  enhanced  damage  in  combined 
environments,  it  would  seem  that  a  surface  preconditioned  by  solid 
particle  impacts  would  provide  numerous  sites  for  the  initiation  of  rain 
erosion  damage.  In  other  words,  the  long  incubation  period  in  Figure  18 


The  rainfield  in  the  AFML/Bell  rotating  arm  facility  was  recalibrated 
and  the  plumbing  system  improved  in  1978.  Deviations  in  the  previously 
stated  and  actual  rainfield  parameters  were  uncovered  at  that  time,  so 
the  erosion  rate  could  be  changed  simply  due  to  the  calibration  differ- 
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(the  duration  of  which  is  dependent  on  the  impact  velocity)  would  be 
by-passed  and  the  steady-state  erosion  rates  for  the  rain  environment 
would  be  rapidly  achieved.  Also  due  to  the  pre-existing  fractures,  the 
effectiveness  of  the  penetration  and  upheaval  mechanism  would  be 
amplified  for  the  rain  drop  collisions.  The  combined  effect  of  exposure 
to  a  dust  and  then  a  rain  environment  should  be  much  more  severe  than 
simply  encountering  a  rain  environment  with  a  pristine  material  surface 
for  a  semi-brittle  material.  On  the  other  hand,  flight  through  rain 
followed  by  exposure  to  dust  should  produce  only  a  slight  increase  in 
the  solid  particle  erosion  rates  due  to  the  localized  nature  of  the 
solid  particle  impact  damage.  At  this  time  these  observations  are 
conjectures  and  require  experimental  verification  and  quantification. 
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SECTION  6 


EXISTING  PARTICLE  IMPACT  TESTING  CAPABILITY 

The  previous  section  provided  an  assessment  of  the  prevailing  data 
base  for  the  particulate  erosion  of  composite  material  constructions. 
The  main  conclusion  from  this  survey  is  that  there  is  essentially  no 
data  available  for  the  types  of  materials  which  are  of  current  interest. 
Another  conclusion  which  is  evident  from  this  survey  is  that  the  number 
of  existing  facilities  which  can  provide  dust  erosion  data  at  impact 
velocities  from  500  to  1500  fps  is  relatively  small  and  the  specimen 
dimensions  are  quite  limited.  The  more  positive  result  from  this  survey 
is  the  clear  indication  provided  in  Section  2  of  how  laboratory-scale 
test  data  can  be  utilized  to  provide  estimates  of  the  extent  of  the 
erosion  damage  on  a  wide  range  of  materials  for  actual  flight  scenarios. 
The  problem  of  how  to  obtain  useful  particulate  erosion  damage 
assessments  will  be  addressed. 

Particulate  erosion  tests  are  usually  designed  to  satisfy  one  or 
more  of  the  following  three  objectives:  straightforward  screening  of 
materials  exposed  to  a  standard  test  condition,  identification  of  damage 
mechanisms  in  order  to  improve  the  erosion  resistance  of  a  particular 
material,  and  the  generation  of  data  in  support  of  a  predictive  modeling 
effort.  The  problem  areas  have  mainly  involved  studying  the  effects  of 
a  single  or  relatively  small  number  of  particle  impacts,  such  as, 
hypervelocity  particle  impacts  on  re-entry  vehicle  nosetips  or  water 
drop  collisions  on  fragile  infrared-transmitting  materials  or  space 
shuttle  tiles. 

Motivated  for  the  most  part  by  the  need  for  evaluating  the 
relative  rate  of  material  removed  due  to  solid  particle  or  water  drop 
impingement  rather  than  in  understanding  the  material  removal  process 
itself,  a  diversity  of  test  configurations  have  been  utilized  to 
generate  the  screening  assessments  required.  Erosion  testing  has 
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generally  been  undertaken  to  satisfy  a  very  specific  need  for  a 
particular  application.  Typically  a  test  method  is  devised 

representative  of  the  operational  environment.  A  material  selection 
procedure  is  pursued  and  then  the  testing  is  terminated.  It  is 
difficult  to  construct  a  coherent  series  of  test  results  on  this  basis. 
On  the  other  hand  there  have  been  a  few  programs  which  have  been  carried 
out  to  investigate  the  erosion  process;  these  have  been  referred  to 
fairly  extensively  in  the  erosion  literature.  Unfortunately,  the 

available  solid  particle  .  erosion  test  results  have  generally  been 
reported  with  insufficient  information  pertaining  to  the  test  conditions 
to  provide  a  valid  comparison  between  the  data  obtained  in  one 
laboratory  with  that  from  another  laboratory.  This  situation  is 
complicated  further  by  the  fact  that  changes  in  the  construction  of  the 
same  test  apparatus  can  significantly  influence  the  delivery  of 
particles  to  the  target  material  as  shown  by  Uuemois  and  Kleis  (1975), 
Wolak ,  et  al.  (  1977)  and  Maji  and  Sheldon  (  1979). 

The  primary  reasons  for  the  discrepancies  in  the  correspondence 
between  the  erosive  response  of  aircraft  components  based  on  erosion 
facility  test  results  and  their  actual  performance  under  operational 
conditions  are: 

•  the  lack  of  adequate  characterization  of  the  particle  (sand  or 
rain)  interaction  with  the  aerodynamic  flow  fields  around  the 
component  during  flight  (recognizing  that  characterization  of 
these  particulate  environments  in  general  is  a  difficult 
problem) , 

•  the  lack  of  characterization  and  verification  of  the  particle 
impact  conditions  within  the  erosion  test  facilities  used  to 
generate  the  erosion  data, 

•  the  lack  of  a  procedure  which  indicates  how  the  erosion  test 
results  relate  to  operational  conditions. 
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The  model  described  in  Section  2  provides  a  general  framework  for 
relating  erosion  test  results  with  the  actual  flight  conditions  and 
erosive  environments  in  order  to  predict  the  material  losses  a 
structural  component  may  experience.  The  model  can  be  applied  to 
material  loss  estimates  for  aircraft  components.  This  model  clearly 
indicates  that  once  accurate  erosion  test  data  is  obtained  it  can  be 
used  to  evaluate  a  broad  range  of  flight  scenarios.  This  is  a  very 
significant  result.  This  approach  affords  the  opportunity  to  use  a 
well-planned,  laboratory-scale  test  program  to  generate  a  data  base  from 
which  an  unlimited  number  of  flight  scenarios  and  trade-off  studies  can 
be  undertaken,  compared  with  simulation  tests  which  have  a  narrow  range 
of  application  and  have  to  be  repeated  as  the  flight  conditions  are 
changed.  This  is  not  to  imply  that  the  suggested  approach  is  easy  to 
implement,  but  once  accomplished  it  has  a  considerably  higher  utility 
and  cost-effectiveness  than  the  alternatives. 

The  need  for  multiple  solid  particle  collisions  at  impact 
velocities  from  500  to  1500  fps  limits  the  facilities  which  can  achieve 
these  particle  impact  requirements  under  reasonably  controlled 
conditions.  The  types  of  facilities  commonly  used  to  obtain  solid 
particle  erosion  data  will  be  described  and  then  the  possibilities  for 
satisfying  the  higher  velocity  requirements  will  be  considered. 

The  most  widely  used  test  configuration  for  solid  particle  erosion 
is  the  sand  blast  apparatus.  Although  some  version  of  this  device  has 
been  used  for  laboratory-scale  erosion  evaluations  for  over  seventy 
years,  it  has  only  been  quite  recent  that  the  test  parameters  associated 
with  particular  types  of  equipment  have  been  investigated.  The  sand 
blast  equipment  used  for  providing  a  multiple  particle  field  ingests 
particulate  matter  into  a  flowing  gas  stream  which  subsequently  impinges 
on  a  test  piece  through  a  small  diameter  tube:  typically  less  than  12  mm 
(0.5  in).  The  impact  velocity  and  angle  of  attack  of  the  particles  in 
the  gas  stream  are  governed  in  part  by  their  size,  shape,  and  density. 
The  maximum  particle  impact  velocities  are  generally  less  than  500  fps. 


A  facility  of  this  type  which  operates  at  higher  particle  impact 
velocities  was  developed  during  the  early  1970s  at  the  Solar  Turbines 
International,  San  Diego,  California.  In  the  Solar  facility  the  dust  is 
injected  into  a  high-pressure  air  stream.  The  length  of  the  delivery 
tube  is  on  the  order  of  10  ft  (3.05  m)  so  the  particles  have  sufficient 
time  in  the  air  flow  to  achieve  their  terminal  velocities.  The 
particles  exit  the  system  through  a  0.375  in  (0.95  cm)  nozzle.  With 
this  arrangement  the  eroded  area  is  0.11  in  (0.71  cm  )  which  may  be 
marginal  with  respect  to  the  construction  of  some  of  the  composites 
which  may  be  evaluated.  The  maximum  air  velocity  is  limited  by  the 
sonic  velocity:  1120  ft/s  (342  ms  ^) .  The  maximum  particle  velocities 
depend  on  the  mean  size  of  the  particles.  The  maximum  particle  velocity 
for  a  10  Um  particle  is  935  ft/s  (285  ms  ^);  the  impact  velocity  in  this 
facility  is  greatly  reduced  for  particles  on  the  order  of  100  (jm.  The 
particle  impact  conditions  again  fall  short  of  the  desired  regime. 

Another  facility  which  provides  higher-velocity  dust  erosion 
impacts  (velocities  up  to  800  fps  for  50  um  particles)  was  recently  put 
in  operation  at  PDA  Engineering,  Inc.,  in  Santa  Ana,  California.  The 
existing  capability  for  this  facility  will  be  briefly  described  as  it  is 
known  at  the  present  time. 

The  particle  impingement  on  the  target  is  accomplished  by 
translating  the  target  across  the  particulate  flow  from  one  to  four 
nozzles.  The  nozzle  inside  diameter  is  0.19  in.  and  the  spacing  between 
the  linear  array  of  nozzles  is  about  0.25  in.  The  nozzle  array  can  also 
be  oscillated  in  an  up  and  down  motion  which  spreads  the  particles 
uniformly  over  a  0.5  in.  high  strip  approximately  1.75  in.  long.  The 
test  specimen  is  mounted  to  a  7  inch  square  plate  that  translates  across 
the  nozzles.  The  motions  of  this  plate  are  computer  controlled,  so  the 
traverse  and  speed  at  which  the  plate  is  moved  can  be  regulated.  Impact 
angles  currently  ranging  from  20  to  90  degrees  can  be  used.  At  impact 
angles  less  than  90  degrees  the  distance  from  the  centerline  of  the 


nozzle  array  to  the  impact  location  on  the  plate  is  maintained  at  a 
constant  value  once  it  is  set. 

The  working  area  for  uniform  particle  flows  is  about  6  inches 
square  in  order  to  account  for  end  effects.  PDA  indicates  that  each 
location  on  the  specimen  receives  the  same  particle  distribution.  The 
edges  of  the  specimen  are  protected  by  placing  a  frame  around  the 
specimen.  The  specimen  mounting  conditions  are  flexible  so  the  edge 
restraints  can  be  designed  to  minimize  the  influence  of  the  specimen 
boundaries  on  the  rate  of  material  removal.  The  maximum  particulate 
exposure  durations  can  range  from  20  to  60  minutes  depending  on  the 
particle  size  and  the  flow  rate  desired. 

The  particles  are  injected  into  the  airflow  from  four  fluidized 
beds  which  are  charged  with  particles  of  a  particular  size.  From  one  to 
four  beds  can  be  operating  depending  on  the  desired  particle  concentra¬ 
tion.  The  particle  concentration  can  also  be  changed  by  changing  the 
flow  rate  in  the  fluidized  bed.  All  four  beds  can  contain  the  same 
particle  sizes  or  the  sizes  can  be  different  in  each  bed  to  provide  a 
prescribed  particle  distribution.  The  particle  velocity  achieved  for 
each  particle  dimension  is  dependent  on  the  particle  size,  so  the  flow 
velocities  have  to  be  adjusted  for  each  bed  if  a  range  of  particle  sizes 
is  to  impact  the  specimen  at  approximately  the  same  velocity  as  would  be 
the  case  for  an  aircraft  flying  into  a  dust  cloud.  The  exactness  of 
this  condition  is  dependent  upon  the  particle  velocities  within  the  dust 
cloud  prior  to  colliding  with  the  aircraft.  The  flow  out  of  each  nozzle 
is  from  a  single  fluidized  bed,  so  mixing  of  particle  sizes  is 
accomplished  by  the  translational  motions  of  the  test  specimen  past  the 
nozzles.  The  particle  velocities  are  determined  using  a  xenon  strobe 
and  photographing  the  particles  along  a  plane  parallel  to  the  flow 
direction. 

Solid  particles  can  also  be  injected  into  a  flowing  gas  stream  as 
in  a  wind  tunnel.  The  specimen  is  then  completely  engulfed  by  the 
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particle-laden  gas  which  offers  the  advantage  of  eliminating  the 
transition  region  around  the  eroded  spot  on  the  test  piece  when  the 
particles  are  propelled  from  a  nozzle  or  tube  as  in  the  sand  blast 
arrangement.  A  system  of  this  type  has  been  described  by  Grant  and 
Tabakoff  (1975)  and  Tabakoff  and  Wakeman  (1979)  at  the  University  of 
Cincinnati.  The  maximum  air  velocity  achievable  in  the  dust  erosion 
tunnel  for  ambient  temperatures  at  the  University  of  Cincinnati  is  450 
ft/s  (140  ms  ■*■) .  Wind  tunnels  have  been  used  in  the  past  but  the 
characterization  of  the  dust  environments  in  the  tunnels  was  a  problem 
in  the  early  investigations.  At  the  present  time  it  does  not  appear 
that  dust-laden  airflows  are  incorporated  into  any  existing  transonic  to 
low  supersonic  wind  tunnel  within  the  United  States. 

Rotating  arm  devices  in  which  the  specimen  is  mounted  at  the  tip 
of  a  rotor  blade  and  then  passes  through  a  stream  of  falling  particles 
are  also  used  for  solid  particle  erosion  evaluations  (Kleis,  1969;  Tilly 
and  Sage,  1970;  Adler,  Morris,  and  Wahl,  1972).  The  rotating  arm 
configuration  has  the  advantage  of  having  the  specimen  impact  a  graded 
particle  distribution  (such  as  specified  by  U.S.  Military  Specification, 
Mil-E-5007C,  which  includes  particle  sizes  over  the  range  from  0  to  1000 
cm  at  a  uniform  velocity.  Such  conditions  are  characteristic  of  high 
speed  structural  components  in  a  dust  field  where  the  particles  are 
traveling  at  considerably  lower  velocities.  The  rotating  arm  configura¬ 
tions  also  have  potential  for  being  operated  at  higher  impact  velocities 
compared  with  sand  blast  and  wind  tunnel  test  facilities,  but  the  sample 
mass  and  dimensions  are  quite  small  and  the  sample  configurations  are 
also  highly  restricted.  Bell  Aerospace  Textron  rotating  arm  can  operate 
at  velocities  of  around  100  ft/s  to  about  3000  ft/s.  However  as  Mach  1 
is  exceeded  the  sample  dimensions  are  limited. 

The  variations  which  exist  in  solid  particle  erosion  test 
apparatus,  even  for  sand  blast  devices,  make  it  difficult  to  discuss 
typical,  or  standard,  test  conditions.  Generally  the  test  results  are 
strongly  related  to  the  test  apparatus  used.  In  addition  the  testing 


procedures  for  solid  particle  erosion  vary  widely  from  laboratory  to 
laboratory,  so  it  is  difficult  to  make  an  absolute  comparison  between 
the  results  obtained  for  one  set  of  test  conditions  to  another.  Also, 
as  more  attention  is  paid  to  the  test  procedures  and  characterization  of 
the  impact  conditions,  additional  factors  which  invluence  the  impacting 
particle/target  interactions  are  being  discovered  (Tilly,  1969;  Tilly 
and  Sage,  1970;  Uuemois  and*  Kleis,  1975;  Grant  and  Tabakoff,  1975; 
Wolak,  et  al.,  1977;  Young  and  Ruff,  1977;  Tabakoff,  Kotwal,  and  Hamed, 
1979;  Maji  and  Sheldon,  1979;  Lapides  and  Levy,  1979).  Therefore  at  the 
present  time  the  test  facility  and  the  test  procedure  must  be  considered 
in  conjunction  with  the  test  results. 


Surveys  of  solid  particle  erosion  testing  facilities  and 
procedures  indicate  there  is  a  definite  need  for  devising  a  facility  to 
provide  an  accurate  measure  of  the  erosive  response  of  materials  and/or 
to  provide  detailed  characterization  of  the  actual  impact  parameters  for 
existing  facilities.  This  is  not  a  simple  undertaking  due  to  the  number 
of  parameters  which  are  associated  with  multiple  collisions  of 
irregular-shaped  particles  on  a  material  surface,  whose  general 
topography  and  properties  are  changing  with  the  length  of  the  exposure. 
The  characterization  process  for  a  single  particle  type  and  material  is 
in  itself  quite  extensive.  Many  factors  have  to  be  evaluated  to  develop 
accurate  measures  of  the  mass  of  particles  actually  contacting  the 
surface,  their  impact  velocity  and  impingement  angle,  and  the  nature  of 
the  particles  delivered  to  the  surface  compared  with  the  initial 
properties  and  size  of  the  particles  (especially  for  soft  particulates). 
Concentration  effects  are  obviously  important.  Particle  concentration 
can  influence  the  extent  of  particle  shielding  at  the  surface  (the 
accumulation  .  of  particles  on  the  surface  which  prevents  the  incoming 
particles  from  making  direct  contact  with  the  target's  surface), 
particle  rebound  characteristics  which  change  as  a  function  of  exposure 
time  also  affect  the  trajectories  of  the  incoming  particles,  and  the 
occurrence  of  particle  embedding  in  the  target's  surface  modifies  the 
measured  erosive  response  of  the  base  material. 


An  overview  of  Che  evaluation  of  the  system  implications  of 
aircraft  flying  ctirougn  nuclear-generated  dust  environments  will  be 
described  in  the  next  section.  The  explicit  erosion  testing 
requirements  will  be  outlined  in  conjunction  with  this  overview. 


SECTION  7 


GENERAL  EXPERIMENTAL  REQUIREMENTS 


The  statement  is  often  made  that  the  erosion  data  obtained  from 
laboratory-scale  erosion  facilities  is  not  applicable  to  the  conditions 
prevailing  lor  actual  flights  through  dusty  environments.  Contrary  to 
this  commonly-accepted  observation,  it  is  our  opinion  that  if  the 
laboratory  data  is  used  in  the  right  way,  very  realistic  predictions  can 
be  made  of  the  erosive  losses  for  a  broad  spectrum  of  aircraft  flight 
scenarios.  This  approach  offers  the  advantage  of  not  having  to  carry 
out  a  significant  number  of  costly  flight  tests  even  if  it  was  possible 
to  undertake  full-scale  testing  in  dusty  environments.  Since  it  is 
highly  impractical  to  fly  aircraft  in  dusty  environments  representative 
of  nuclear-generated  dust,  the  suggested  laboratory-scale  testing 
provides  a  viable  solution.  The  accuracy  of  the  erosive  loss 
predictions  is  merely  dependent  on  the  level  of  effort  devoted  to  the 
adjunct  analyses  required  to  transform  the  laboratory  erosion  data  to 
the  aircraft  flight  conditions.  The  statements  made  in  Section  2  are 
reiterated,  since  it  should  be  understood  from  the  outset  that  once 
reliable,  laboratory  erosion  data  is  obtained  for  the  range  of 
conditions  which  will  be  specified,  this  data  remains  valid  for  all 
flight  scenarios  and  environmental  conditions  which  may  be  considered. 

The  analyses  and  observations  from  Sections  2  and  6  provide  the 
background  for  a  meaningtul  approach  for  an  assessment  or  the  potential 
dust  erosion  damage  to  aircraft  flying  at  high  subsonic  to  low 
supersonic  speeds.  An  overview  or  how  such  an  assessment  can  be 
accomplished  rs  outlined  in  Figures  19  and  20.  The  interrelationships 
in  Figure  19  provide  a  clear  perspective  tor  tue  main  features  of  the 
dust  erosion  assessment.  The  output  from  the  analyses  and  experimental 
data  generation  will  be  a  computational  model  for  describing  the  extent 
of  the  material  removal  which  may  occur  for  an  arbitrary  flight  plan. 
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Figure  20.  Prediction  of  dust  erosion  effects  on  aircraft  components 
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The  resulting  computational  model  may  be  used  to  examine  a  wide  range  of 
flight  conditions. 


The  primary  constituents  of  ttie  computational  model  listed  in 
Figure  19  will  be  described.  It  is  important  to  note  that  estimates  of 
the  mass  removal  can  be  made  using  a  minimal  amount  of  information.  As 
the  inputs  to  the  general  approach  are  refined  the  uncertainties  in  the 
erosion  predictions  will  correspondingly  be  reduced. 

Either  a  superficial  or  detailed  description  of  the  erosive 
environment  can  be  used  to  bracket  the  particulate  impact  conditions 
wnich  w.ll  be  significant.  The  initial  survey  of  the  particle  impact 
parameters  is  used  to  limit  the  number  of  particle  impact  conditions 
which  tiave  to  be  considered.  However  subsequent  modification  of  the 
dust  environments  will  not  alter  the  erosion  test  conditions.  They 
should  oe  selected  to  be  broad  enough  to  accommodate  the  anticipated 
dust  environments.  Since  a  very  detailed  and  complete  erosion  test 
matrix  can  become  overwhelming  from  ttie  standpoint  of  the  scope  of  the 
test  program  required,  judgment  has  to  be  exercised  as  to  the  conditions 
which  will  provide  a  reasonable  test  plan  as  well  as  a  reliable  data 
base.  With  proper  selection  of  the  test  conditions  the  data  base  can  be 
expanded  with  extrapolations  of  the  experimental  results. 

The  motivation  for  elaborating  the  utility  of  the  experimental 
erosion  data  base  is  that  it  is  sometimes  thought  that  the  high  level  of 
uncertainties  associated  with  specifying  the  nuclear-generated  dust 
environments  should  be  reduced  before  any  effort  be  devoted  to 
developing  a  much  more  representative  and  reliable  erosion  data  base 
than  exists  at  the  present  time.  While  significant  discrepancies  in  the 
impact  parameters  derived  from  the  lofted  dust  and  debris  models  and  tne 
dispersion  of  this  dust  in  ttie  atmosphere  will  alter  the  predictions  of 
aircraft  erosion,  an  initially  conservative  estimate  can  be  provided  by 
assuming  a  worst-case  condition.  If  the  resulting  levels  of  material 
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removal  are  not  significant,  the  investigation  of  the  erosive  effects 
does  not  have  to  be  pursued  further.  If  the  resulting  predictions  of 
material  loss  are  marginal  or  unacceptable,  further  refinement  of  the 
dust  environments  to  be  encountered  would  be  warranted.  As  the 
nuclear-generated  dust  models  are  refined,  specification  of  the  revised 
dust  environments  can  be  readily  and  easily  incorporated  into  the 
erosion  assessment.  A  much  more  important  consideration  for 
implementing  the  suggested  approach  is  to  provide  reliable  laboratory 
erosion  data. 

In  order  to  obtain  reliable  erosion  data,  the  particle  impact 
conditions  within  the  laboratory  erosion  facility  and  the  particle 
impact  conditions  at  the  surface  of  the  vulnerable  aircraft  components 
should  be  well-character ized .  This  aspect  of  utilizing  the  laboratory 
erosion  data  is  usually  neglected  and  leads  to  the  general  lack  oi 
correspondence  between  the  laboratory  results  and  the  actual  flight 
conditions.  In  the  approach  shown  in  Figures  19  and  20  it  is  very 
important  that  the  aerodynamics  of  incoming  particles  impacting  t lie  test 
specimens  in  the  erosion  test  facility  and  for  the  aircraft  be  defined. 
In  the  erosion  facility  this  can  be  accomplished  using  diagnostic 


equipment  to  directly  observe  the  particle  impacts  in  addition  to 
analytically  defining  their  trajectories.  For  the  aircraft  the  particle 
trajectories  will  more  likely  ha-ve  to  be  specified  from  analytical  or 
computational  models.  "The  difference  in  these  particle  impact 
conditions  has  been  explicitly  introduced  in  Equation  (5).  The  level  to 
which  these  particle  impact  conditions  can  be  prescribed  will  determine 
the  accuracy  of  the  predictions  which  can  be  made  by  this  approach.  The 
general  relations  for  incorporating  the  experimental  data  and  for 
predicting  the  mass  removal  for  actual  flight  scenarios  involves  the 
explicit  evaluation  of  in  Equations  (5)  and  (12)  in  Section  2. 

Since  a  complete  description  of  the  mass  loss  as  a  function  of  the 
particle  type,  size,  impact  angle,  impact  velocity,  and  particle 
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concentration  tor  a  single  material  involves  hundreds  or  particle  impact 


i 

i 

1 
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tests,  discretion  has  to  be  exercised  in  the  number  of  impact  tests 
which  are  reasonable.  If  a  number  of  materials  are  to  be  evaluated,  a 
baseline  material  can  be  selected  and  fairly  extensively  characterized. 
The  test  program  for  these  additional  materials  can  be  limited  to  the 
extent  that  is  necessary  to  at  least  establish  the  trends  (with  respect 
to  impact  angle,  particle  size  and  impact  velocity)  which  are  adequate 
to  define  their  particulate  erosive  behavior  relative  to  the  baseline 
material.  A  limited  number  of  test  conditions  can  be  prescribed  for  the 
new  materials  of  the  same  general  type  in  the  test  matrix  and  these  data 
points  superposed  on  the  more  complete  curves  for  the  baseline  material. 
If  the  same  general  trends  appear  as  for  the  baseline  material,  a 
limited  number  of  tests  will  be  sufficient  to  compare  each  new  material 
with  the  baseline.  If  anomalies  are  noted  then  more  tests  (a  second 
interaction)  will  be  required  to  describe  the  behavior  of  the  new 
material. 

In  addition  to  obtaining  quantitative  measures  of  the  material 
removal  rates  the  test  program  also  offers  the  advantage  of  providing 
specimens  tor  materials  characterization.  Microscopic  examination  of 
the  damage  produced  from  the  onset  of  material  removal  to  accelerated 
rates  of  material  removal  provides  significant  information  on  how  the 
erosion  resistance  of  a  particular  material  can  be  improved.  Identi¬ 
fication  of  the  erosion  mechanisms  can  be  accomplished  by  removing 
specimens  at  various  stages  of  the  erosion  process.  These  can  then  be 
used  to  construct  a  physical  description  of  the  microscopic  mechanisms 
responsible  for  the  onset  and  subsequent  progression  of  damage. 

In  order  to  gain  some  appreciation  for  the  magnitude  of  the 
erosion  test  program  a  representative  test  plan  will  be  outlined. 
Suppose  only  one  type  of  dust  is  considered,  for  example,  crushed 
quartz.  The  silica  can  be  sieved  in  accordance  with  the  standard  Tyler 
sieve  openings  beginning  witii  38  urn.  Four  size  ranges  for  the  particle 
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sizes  are  reasonable  spanning  Che  dimensions  from  38  co  180  urn.  Since 
che  impact  angle  dependence  of  the  material  removal  rates  for  the 
candidate  composites  is  not  known,  the  impact  angles  of  15,  25,  45,  60, 
and  90  degrees  would  have  to  be  included  for  the  baseline  material  at 
least.  The  evaluation  of  the  impact  velocity  dependence  of  the  material 
removal  rates  would  require  a  minimum  of  three  impact  velocities.  There 
are  now  sixty  particle  impact  conditions  specified  for  determining  the 
erosion  coefficients,  u  in  Equation  (5),  for  the  baseline  material. 
However  the  evaluation  of  the  steady-state  erosion  rates  requires  a 
series  of  repetitive  exposures  for  each  particle  impact  condition  in 
order  co  obtain  the  curve  illustrated  in  Figure  1 .  Possibly  eight  to 
nine  exposures  are  required  in  the  erosion  facility  for  this  purpose. 
The  specimen  has  to  be  accurately  weighed  after  each  exposure  and 
correlated  with  the  mass  of  particulates  impacted  during  each  successive 
exposure  increment.  Thus  for  the  baseline  material  possibly  500 
incremental  exposures  to  the  dust  environment  are  necessary  to  fully 
evaluate  the  -erosion  coefficients.  At  least  triplicate  specimens 
exposed  to  the  same  particle  impact  conditions  are  required  to  provide 
adequate  confidence  .in  the  resulting  erosion  coefficients.  How 
expediently  this  could  be  accomplished  depends  on  the  configuration  of 
the  facility  used  for  the  dust  erosion  testing. 


Once  the  erosion  coefficients  have  been  evaluated,  the  particle 

impact  conditions  have  been  properly  characterized,  and  the  erosion 

model  has  been  programmed  for  digital  computations  predictions  of  the 

erosion  of  aircraft  materials  can  be  made.  This  is  the  scheme  shown  in 

Figure  19.  The  erosion  model  is  tnen  combined  with  a  definition  of  the 

dust  environment  which  can  be  used  as  input  data  for  the  model  and 

evaluations  of  the  aerodynamic  flow  fields  over  the  aircraft  component 

* 

of  interest  as  indicated  in  Figure  20.  Predictions  of  the  extent  of  the 
material  removal  from  this  component  can  then  be  provided. 
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Without  the  option  of  instrumenting  an  aircraft  and  flying  through 
a  well-characterized  dust  environment,  a  more  feasible  assessment  or  the 
accuracy  of  the  erosion  predictions  for  the  aircraft  component  may  be 
obtained  with  wind  tunnel  tests  for  at  least  a  section  of  the  component. 
The  general  plan  for  conducting  the  wind  tunnel  tests  is  illustrated  in 
Figure  21.  The  erosion  model  would  be  evaluated  for  the  conditions 
within  the  wind  tunnel.  The  predictions  from  the  model  would  then  be 
compared  with  the  post-test  measurements  of  the  material  removal  from 
the  actual  aircraft  component  after  a  sufficiently  long  exposure  to  the 
particle  impact  conditions  within  the  wind  tunnel. 

While  the  effort  described  in  this  section  is  substantial,  it 
provides  a  rational  approach  for  evaluating  the  dust  erosion  effects 
which  in  the  long  run  will  be  cost-effective  and  will  apply  to  a  wide 
range  of  flight  conditions. 


AIRCRAFT 

COMPONENT 


Figure  21.  Comparison  of  erosion  model  predictions  with 

simulated  aircraft  component  flight  conditions. 
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CONCLUSIONS 


The  existing  solid  particle  erosion  data  for  non-metallic 
composite  materials  has  been  summarized  and  reviewed.  This  survey 
indicates  that  the  erosion  data  for  the  particle  sizes  and  impact 
velocities  of  interest  is  essentially  not  available.  The  limited  number 


of  erosion  rates  which  are  available  have  been  used  as  the  basis  for  a 


very  preliminary  estimate  of  the  material  removal  which  may  be 
experienced  by  aircraft  flying  through  dusty  environments.  This  is 
certainly  a  highly  speculative  approach,  but  it  represents  the  best 
estimates  chat  can  be  made  at  the  present  time.  The  available  data 
contains  hign  levels  of  uncertainty,  since  it  is  for  very  ill-defined 
composite  materials  and  the  test  procedures  used  may  introduce 


signiticant  errors. 


An  erosion  model  which  clearly  indicates  the  erosion  and  environ¬ 
mental  parameters  required  and  how  they  affect  the  material  removal 
rates  was  used  to  predict  the  magnitude  of  the  material  removal  which 
may  occur  for  hypothetical  flight  scenarios.  The  indication  from  this 
initial  effort  is  that  the  solid  particle  material  removal  rates  for 
advanced  composites  can  be  significant. 


This  investigation,  which  represents  the  best  material  loss 


estimates  that  can  be  made  currently,  indicates  that  a  properly  directed 
experimental  effort  is  required  to  improve  the  reliability  of  the  damage 
assessment  of  composite  aircraft  components  exposed  to  particulate 


environments. 
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DEPARTMENT  OF  ENERGY 

SANDIA  NATIONAL  LABORATORIES 

ATTN:  TECH  LIB  3141  RPTS  RCVG  CLRK 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  OSWR/NED 

U  S  ARMS  CONTROL  &  DISARMAMENT  AGCY 
ATTN:  H COOPER 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

BOEING  CO 

ATTN:  ADACOSTA 

BOEING  MILITARY  AIRPLANE  CO 
ATTN:  DSAWDY 

CALSPAN  CORP 

ATTN:  C  PADOVA 
ATTN:  M  DUNN 

CARPENTER  RESEARCH  CORP 
ATTN:  HJ  CARPENTER 

DAMASKOS,  INC 

ATTN:  N DAMASKOS 

DAYTON,  UNIVERSITY  OF 
ATTN:  B  WILT 


GENERAL  RESEARCH  CORP 
ATTN:  DMIHORA 
ATTN:  R CRAWFORD 
6  CYS  ATTN:  W  ADLER 

IRT  CORP 

ATTN:  D  WOODALL 

KAMAN  SCIENCES  CORP 
ATTN:  D COYNE 
ATTN:  LMENTE 
ATTN:  R  RUETENIK 
ATTN:  W  LEE 

KAMAN  SCIENCES  CORP 
ATTN:  E  CONRAD 

KAMAN  TEMPO 

ATTN:  DASIAC 

KAMAN  TEMPO 

ATTN:  DASIAC 

LOCKHEED  CORPORATION 
ATTN:  ASCHUETZ 
ATTN:  B OSBORNE 

LOCKHEED-GEORGIA  COMPANY 
ATTN:  ABONDERUD 
ATTN:  R  KELLY 
ATTN:  T  HODGSON 

MCDONNELL  DOUGLAS  CORP 
ATTN:  H  SAMS 
ATTN:  L  KOCH 
ATTN:  R  J  PANDORF 

MCDONNELL  DOUGLAS  CORP 
ATTN:  J  TRACY 
ATTN:  L  COHEN 

MCDONNELL  DOUGLAS  CORP 
ATTN:  J  GAIDULIS 
ATTN:  JMCGREW 
ATTN:  M  POTTER 

NORTHROP  CORP 

ATTN:  BAGARWAL 
ATTN:  G KENDALL 

PACIFIC-SIERRA  RESEARCH  CORP 

ATTN:  H  BRODE,  CHAIRMAN  SAGE 

R  &  D  ASSOCIATES 

ATTN:  AKUHL 
ATTN:  FA  FIELD 
ATTN:  P RAUSCH 

RAND  CORP 

ATTN:  P  DAVIS 


DEPT  OF  DEFENSE  CONTRACTORS  (CONTINUED)  SCIENCE  APPLICATIONS  INTL  CORP 

ATTN:  AMARTELLUCCI 

RAND  CORP 

ATTN:  B  BENNETT  SOUTHERN  RESEARCH  INSTITUTE 

ATTN*  C  PEARS 

ROCKWELL  INTERNATIONAL  CORP  ATTN:  S  CAUSEY 

ATTN:  A  MUSICMAN 

ATTN:  P  MASON  TOYON  RESEARCH  CORP 

ATTN:  SMELLINS  ATTN:  B  GRAGG 

ATTN:  J  CUNNINGHAM 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J  COCKAYNE 


